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Abstraci-VJe  describe  a  new  anomaly  detection  algorithm  based 
on  Electrical  Impedance  Tomography  (EIT)  technique.  When 
only  the  boundary  current  and  voltage  measurements  are 
available,  it  is  not  practically  feasible  to  reconstruct  accurate 
high-resolution  cross-sectional  resistivity  images  of  a  subject.  In 
this  paper,  we  focus  our  attention  on  the  detection  of  the  location 
and  size  of  anomalies  with  resistivity  values  different  from  the 
background  tissues.  We  show  the  performance  of  the  algorithm 
from  experimental  results  using  32-channel  EIT  system  and 
saline  phantoms.  The  algorithm  is  applicable  to  the  detection  of 
cancerous  tissues  in  the  breast. 
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I.  INTRODUCTION 

Electrical  impedance  tomography  (EIT)  is  a  medical 
imaging  technique  to  produce  cross-sectional  resistivity  (or 
conductivity)  images  of  the  human  body  [1],  [2].  Electrodes 
are  placed  on  the  body  surface  to  inject  currents  and  measure 
voltages.  Several  image  reconstruction  algorithms  are 
published  to  image  the  spatial  distribution  of  electrical 
resistivity  [2].  Since  1980s,  this  technique  has  been  received 
considerable  attention  in  medical  imaging  area  because  it  may 
visualize  totally  new  physiological  information  with  many 
possible  clinical  applications.  In  particular,  considerable 
efforts  have  been  devoted  to  applications  in  cardiac, 
respiratory,  and  breast  imaging. 

In  EIT,  we  need  to  find  the  reconstruction  map  from  the 
measured  boundary  data  to  the  cross-sectional  resistivity 
distribution.  Unfortunately,  unlike  the  other  medical  imaging 
modalities  such  as  xray  CT  and  MRI,  EIT  suffers  from  the 
nonlinearity  between  the  measured  boundary  data  and  the 
internal  resistivity  distribution.  Moreover,  reconstruction 
algorithms  are  severely  ill-posed  due  to  the  inherent  low 
sensitivity  of  boundary  measurements  to  any  changes  of 
internal  tissue  resistivity  values.  This  ill-posedness  and 
nonlinearity  make  EIT  not  so  successful  in  producing  high- 
resolution  images.  For  the  present,  with  a  limited  set  of 
measured  boundary  data,  imaging  high-quality  resistivity 
distribution  seems  to  be  extremely  difficult. 

In  this  paper,  we  focus  our  attention  on  the  detection 
instead  of  imaging.  We  believe  it  is  more  practical  to  extract 
only  necessary  core  information  by  utilizing  all  available 
information  on  the  resistivity  distribution  instead  of  trying  to 
visualize  precise  resistivity  image  extravagantly.  For 
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electrical  impedance  imaging  of  breast  cancer  as  an  example, 
the  resistivity  of  cancerous  tissues  differs  from  that  of 
surrounding  normal  tissues.  In  this  case,  the  size  and  location 
of  the  cancerous  tissues  will  be  the  core  information  to  search 
for  rather  than  trying  to  figure  out  accurate  shape  of  them. 

We  will  formulate  a  model  problem  of  detecting  size  and 
location  of  anomalies  such  as  breast  cancer  by  electrical 
impedance  imaging  technique.  Our  detection  algorithm  is 
based  on  the  recent  mathematical  results  that  introduced  real¬ 
time  three-dimensional  algorithm  to  detect  size  and  location 
of  anomalies  [3],  [4].  Here  we  assume  that  the  differences 
between  resistivity  values  of  anomalies  and  background  are 
known  to  be  in  a  fixed  range.  In  this  method,  the  location  of 
anomaly  is  immediately  determined  by  the  pattern  of  a  simple 
weighted  combination  of  injection  currents  and  measured 
voltages.  We  perform  several  experiments  on  a  saline 
phantom  to  demonstrate  the  accuracy  of  our  detection 
algorithm. 

The  experimental  results  show  that  our  algorithm  is  quite 
accurate  and  stable  against  measurement  noise.  The  detection 
process  is  fast  enough  to  carry  out  all  calculations  using  an 
ordinary  modern  personal  computer  for  real-time  monitoring. 

II.  METHODOLOGY 
A.  Total  Size  Detection 

We  are  interested  in  the  inverse  problem  of  finding 
multiple  anomalies,  \jm._xDj  located  in  a  given  conducting 

body  («=2,  3)  whose  conductivity  distribution  is  of 

the  form, 

<*(*)=i+£wcd,  (1) 

7=1 

where  ^  is  the  characteristic  function  of  Dj  and  q  is  a 

positive  constant  meaning  the  ratio  of  conductivity  values 
between  the  anomalies  and  background.  Then,  for  a  given 

injection  current  g  e  ^(30)  :=|(p  e  L2(d£l) :  JNp  =  o]  ,  the 

electric  potential  u  is  governed  by  the  following  Neumann 
problem: 


V-(aVt/)  =  0  in  Q 

dn  f 

—  =g  on  dQ.  and  ud'k  =  0 

dv  Jan 


(2) 
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where  v  is  the  outward  unit  normal  vector  on  the  boundary 
and  <ZZ  denotes  the  surface  measure. 

In  many  biomedical  applications,  it  may  be  very  desirable 
to  find  the  total  size  using  a  stable  procedure.  One  important 
example  is  the  breast  cancer  detection  problem.  Now,  we 
describe  an  algorithm  to  find  the  total  size  of  multiple 
anomalies  v*J=lDj  •  Suppose  that  uMj=lDj  c=  BR  ,  which  is  the 

ball  with  radius  R  centered  at  x  =  0  contained  in  O.  We 
choose  an  injection  current  g-a-\  where  a  is  a  unit 

constant  vector.  There  exists  a  unique  r ^  0  <  r0  <  R ,  so  that 

U"-v>^=°  (3) 

where  u  and  vr  are  solutions  of  the  Neumann  problems 
9T[Z),g]  and  ,  respectively.  Then,  the  total  size  of 

anomalies  \^jM  d  is  very  close  to  the  size  of  the  ball  b  • 

7-1  7  ro 

Various  numerical  experiments  indicate  that  our  algorithm 
gives  a  nearly  exact  estimate  for  arbitrary  multiple  anomalies 
with  quite  general  conductivity  values.  Its  theoretical 
justification  has  been  studied  in  [3]  where  the  size  estimation 
has  been  proved  for  the  case  where  D  is  a  union  of  small 
disks  Dj,j=  1, . . . ,  M  with  the  following  properties; 

Dj  =Br  (  Zj'j,  r.  <  8,  dist(Z)  ,DJ )  >  L  for  i  =£  j,  dist(Z)y  >  L  (4) 

where  L  is  a  fixed  positive  number  and  8  is  a  small  number. 
Then,  for  the  unique  rQe  (0,7?)  which  minimizes  the  energy 
[3],  the  total  size  of  d  =  can  be  approximated  by 

M 

IN=W+°<e4>-  (5) 

7=1 

B.  Location  Detection 

The  location  search  algorithm  is  based  on  simple  aspects  of 
the  function  Z/(;g,/)  in  Eq.  (6),  which  is  computed  directly 

from  the  boundary  data,  that  is  the  applied  current  g  and  the 
measured  voltage  f. 


Integrating  by  parts  and  using  the  transmission  boundary 
condition  on  dD  [5],  we  get 


77(x;g,/)  +  |aJoV,,<l>(x  -y)  -Vu(y)dy 


w(x)forxe£2  ^ 
OforxelTAft. 


Our  problem  s  now  reduced  to  determine  D  from  the 
known  function  H(x;gf).  How  to  extract  the  location 
information  from  the  representation  formula  in  Eq.  (8)?  In  the 
recent  paper  [4],  we  presented  the  following  asymptotic 
formula  for  the  current  flow. 

Let  £2  be  a  C2  domain  in  R”  and  g^a-v  be  given  for  a 
unit  vector  a  as  a  Neumann  data  for  Eq.  (2).  Suppose 
D  =  u^ZT  satisfies  that  Dj=Br(zj)  with  Tj  <8  , 

dist (Dj  ,Dk)  >  L  for  j^k,  and  dist(Dy.  ,d£l)>  L  for  all  j=  1, 

...,  M  where  8  and  L  are  fixed  positive  numbers.  Then,  the 
gradient  vector  field  of  the  solution  u  in  Eq.  (2)  on  D  is 
approximated  by 


Vw(x)  = 


77  +  JLl 


-  +  0(8  ),  XG  D 


(9) 


By  substituting  the  relation  in  Eq.  (9)  to  Eq.  (8),  the 
formulation  in  Eq.  (8)  is  reduced  to 


H(x\g,f)  = 


It  r  (x-y)-a 
G)>  +  p)J°  | x-y\n 


dy  +  Q(  84) 


(10) 


where  co^  denotes  the  volume  of  a  unit  ball  in  R"  and 

xeRn\h. 

We  can  find  a  key  information  on  the  location  of  D  from 
observations  of  the  function  H(;g,f)  on  R”\3Q  .  If  we 

assume  that  Ci=Bx  c=  R2  and  D  :=Br(z)  is  a  small  ball 

centered  at  z  with  radius  r  in  £2,  the  mean  value  theorem  for 
harmonic  functions  applied  to  Eq.  (10)  yields 


H(x;gJ) 


JLX r2  ( x-z)-a 
2  +  fi  x-z2 


+  <9(s4),xeR2\a 


(11) 


:=L^X(/)/(y)^  ~L ^-y)s(y)dsy  (6) 


where  xeR”\3f2  and  d>  is  the  fundamental  solution  of 
Laplacian  given  by 


3>(x-y)H 


-Llog\X-y\ 
-1  1 


for  n-  2 


4jt  \x-y\ 


(7) 


Applying  the  injection  current  g=(0,l)-v  and  taking  two 
observation  regions  of  =|x=(2,x2)e  K2  \Q]  and 
E2=|x  =(x1,-2)eK2\a}  ,  we  can  derive  the  following 
equation, 


f  .  ,2 

|Llr  x2  -  z2 

_  2+  (J.  (2-  z)2  +  (x2-z2)2 
|Ll  r2  -2  -  z2 

2+  (I  (jq  -Zj)2  +(2  + z2)2 


if  xe  Z, 
if  xg  Z2. 


(12) 


for  77  —  3. 


Therefore,  by  taking  a  zero  point  of  H(x;g,f)  on  £  and  a 
critical  point  of  H(x;g,f )  °n  Z2  ,  we  can  determine  the 
location  of  D.  Fig.  1  illustrates  the  relation  between  the 
location  of  D  and  the  observation  function  H(x\gJ)  on  Z, 
and  Z2 . 


Fig.  1.  Location  of  an  anomaly  and  the  pattern  of 
H(x  ;g ,/ )  for  |i  >  0 


III.  RESULTS 

A.  EIT  Hardware  and  Phantom 

We  present  experimental  results  in  detecting  the  location 
and  size  of  anomalies  in  a  circular  saline  phantom  with  32 
electrodes  attached  on  the  boundary  with  equal  distance 
between  them.  Fig.  2  shows  the  phantom  and  32-channel  EIT 
system.  The  EIT  system  provides  a  way  to  inject  different 
patterns  of  currents  and  measure  voltages  on  the  electrodes.  It, 
therefore,  includes  32  independent  current  sources  whose 
output  is  50  kHz  sinusoid  with  variable  amplitude.  The  EIT 
system  is  interfaced  to  a  PC  and  a  control  program  sets  the 
amount  of  injection  current  for  each  current  source. 


Fig.  2.  EIT  hardware  and  phantom 


While  we  inject  a  pattern  of  currents  through  all  32 
electrodes,  boundary  voltage  on  each  electrode  is  sequentially 
measured.  Voltage  measuring  circuits  consist  of  a  narrow- 
band  variable -gain  ac  amplifier,  a  phase-sensitive 
demodulator,  and  a  12-bit  A/D  converter. 

We  fill  the  phantom  with  electrolyte  solution  of  known 
conductivity  value.  Anomalies  with  different  conductivity 
values  and  shapes  are  placed  inside  the  phantom.  For 
numerical  experiments,  the  circular  phantom  is  normalized  to 
a  unit  disc  £2  in  M2  and  apply  the  injection  current  pattern 
(Neumann  data);  g  =(0,1)  -v(r.) where  r  denotes  electrodes, 

y  =  0,  ...,31. 

B.  Experimental  Results 

Unavoidable  noise  comes  from  various  sources;  the  circular 
phantom  is  not  exactly  a  disc,  the  electrodes  are  not  exactly 
equally  spaced,  and  there  exist  errors  in  the  amount  of 
injection  currents  and  measured  voltage  data.  In  order  to  find 
the  accuracy  of  our  EIT  system,  we  define  the  relative  error 
as 

E ' 

where  /  is  the  analytically  computed  boundary  voltage,  fm  is 
the  measured  boundary  voltage,  andg  is  the  injection  current. 
We  set  the  magnitude  of  injection  current  at  each  electrode  as 
sin(/0)  where  j  is  the  electrode  number  from  0  to  3 1  and  0  is 
the  angle  between  adjacent  two  electrodes.  For  this  kind  of 
injection  current  pattern,  measured  voltages  on  each  electrode 
should  have  the  same  sin(/0)  pattern  and  the  relative  error 
was  4.6%. 


Fig.  3.  Experimental  result  of  the  total  size  detection 


We  put  four  anomalies  in  the  phantom  with  290mm 
diameter  as  shown  in  Fig.  3.  Using  boundary  voltage 
measurements,  the  total  size  of  four  anomalies  (diameter  of 
the  equivalent  ball)  was  computed  as  13 1 .6mm.  Since  the  true 


(b)  Boundary  Measurements 


total  size  is  125.0mm,  the  error  in  size  detection  is  about 
5.3%. 

In  Fig.  4,  we  put  one  anomaly  in  the  phantom  Using  the 
boundary  voltage  measurements  for  the  injection  current 
pattern  of  sin(/'0),  we  found  the  location  of  the  anomaly  as 
shown  in  Fig.  4(c). 


(a)  Target  (b)  Boundary  Measurements 


(c)  Location  detection 


Fig.  4.  Experimental  result  of  the  location  detection 


V.  Conclusion 

We  developed  a  new  anomaly  detection  algorithm  using  an 
EIT  system.  Assuming  that  we  measure  the  boundary  current 
and  voltage  data  on  the  surface  electrodes,  the  detection 
algorithm  finds  the  location  of  an  anomaly  when  the 
resistivity  of  the  anomaly  is  different  from  that  of  the 
background.  It  also  produces  the  total  area  of  multiple 
anomalies. 

When  there  are  multiple  anomalies,  it  finds  the  location  of 
the  anomalies  in  terms  of  the  center  of  resistivity  values.  We 
plan  to  develop  a  new  EIT  system  with  an  improved  accuracy 
and  an  electrode  interface  unit  for  applying  electrodes  on  the 
breast.  This  kind  of  breast  cancer  detection  system  will  be  a 
good  alternative  to  x-ray  mammography  for  patient  screening 
purpose. 
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IV.  DISCUSSION 

Since  the  detection  accuracy  is  dependent  on  the 
measurement  accuracy,  we  need  to  improve  the  performance 
of  our  EIT  system.  We  found  that  the  accuracy  in  the  detected 
size  is  not  dependent  on  the  location  of  anomaly.  This  means 
that  the  detection  of  a  small  anomaly  at  the  center  is  not 
particularly  more  difficult  than  the  detection  of  the  same 
anomaly  at  the  periphery.  When  there  are  multiple  anomalies, 
current  algorithm  finds  the  total  size  and  the  location  in  terms 
of  the  center  of  conductivity  values. 

Future  study  should  include  the  individual  detection  of 
multiple  anomalies.  We  may  use  the  algorithm  described  in 
this  paper  to  find  the  total  size  and  the  location.  Then,  we  can 
use  this  information  as  an  initial  guess  for  an  iterative 
resistivity  image  reconstruction  algorithm.  During  image 
reconstruction,  we  can  also  use  this  information  as  a 
constraint  on  the  reconstructed  resistivity  image. 


